Condensation of phonons in an ultracold Bose gas 
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We consider the generation of longitudinal phonons in an elongated Bose-condensed gas at zero 
temperature due to parametric resonance as a result of the modulation of the transverse trap fre- 
quency. The nonlinear temporal evolution with account of the phonon-phonon interaction leads 
self-consistently to the formation of the stationary state with the macroscopic occupation of a single 
phonon quantum state. 



PACS numbers: 03.75.Kk, 03.75.Nt, 67.40.Db 

The problem of the Bose-Einstein condensation of exci- 
tations is one of most interesting aspects of the condensed 
matter physics. It is well known that under condition of 
thermal equilibrium the chemical potential of excitations 
vanishes and, as a result, their condensate does not form. 
The only way to overcome this obstacle is a creation of 
the stationary state with the conserved number of exci- 
tations in a non-equilibrium system. In the present pa- 
per we describe a mechanism leading to the formation of 
such state for phonon excitations in the Bose-condensed 
atomic gases. The mechanism is based on generating 
phonons with the help of parametric resonance. In the 
Bose-condensed gases the parametric resonance (PR) is 
responsible for the effective energy transfer between dif- 
ferent branches of excitations and temporal evolution of 
the system. In particular, the damping of transverse os- 
cillations in a elongated cylindrical trap at zero temper- 
ature is due to PR with the production of longitudinal 
phonons [jj. The nontrivial picture of the damping in 
these conditions is revealed in the remarkable work of 
the Paris group . The specific features of this picture 
is a result of the nonlinear evolution accompanying PR 
in such system 0. 

We analyze the temporal evolution of interacting lon- 
gitudinal (sound) phonons in an elongated cylindrical 
parabolic trap subjected to permanent modulation of the 
transverse trap frequency loj_. This modulation brings 
about the oscillations of the condensate density and, ac- 
cordingly, sound velocity. Due to PR this results in the 
generation of pairs of sound phonons with the opposite 
momenta and energies close to a half of the modula- 
tion frequency loq/2. In fact, phonon pairs are produced 
within the finite energy interval loq/2 ± Eq. The pa- 
rameter Eo is connected with the modulation amplitude 
I&j^/ojjJ <C 1 by the relation Eq ~ luq ■ \6lj±/iiJ±\. If Eq 
exceeds the damping factor 7 for longitudinal phonons, 
their occupation numbers begin to increase exponentially 
in time. In this case the phonon-phonon interactions 
become significant. As a result, the evolution scenario 
changes drastically. Let us assume that, due to finite 
value of longitudinal size L, a single phonon level alone 
lies within the narrow energy interval 2Eq. In this case 
the main effect of the phonon-phonon interactions is an 



effective renormalization of the level. As is shown below, 
in the course of the self-consistent nonlinear evolution 
the renormalized level reaches the left or right edge of 
the parametric resonance interval in which the growth of 
the phonon occupation numbers ceases. As a result, the 
stationary state of phonons with the macroscopic occu- 
pation numbers in a single quantum state takes place. In 
the new state the steady periodic space modulation of the 
longitudinal density of a Bose gas occurs with the ampli- 
tude proportional to the number of condensed phonons. 

Let us consider the Bose-condensed gas at T = in a 
trap of the cylindrical symmetry with L 3> R, where R is 
the radius of the condensate. Neglecting the edge effects, 
we can write the general equation for the field operator 
^(r, z, t) of atoms in the form 
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where Uq = Anh 2 a/m, a being the s-scattering length. 
Note, considering a dilute and cold Bose gas with domi- 
nant binary collisions we imply as usually that the radius 
of interaction region tq <C a, the gas parameter na 3 <C 1, 
and the correlation length £ > n -1 / 3 , n is the atomic 
density. In this case, the replacement of two-body in- 
teraction with the effective contact potential is relevant. 
The trap transverse frequency depends on time as 

u)(t) = wx(l + rjsinwot), (2) 

According to the results , we introduce the spatial scal- 
ing parameter b(t) and new variables ~~p = r/b(t), r(i). 
The field operator in terms of new variables can be writ- 
ten as 

x(~P,z,t) ft = mr 2 b 

' b' 



* = 



b(t) ' 2h b' (3) 

Substituting these expressions into Eq.Q), we obtain 
straightforwardly 
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provided that functions b(t) and r(t) satisfy the equations 
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Treating the condensate wave function as independent of 
z, we find that Eq.Q expressed in variables p, t describes 
the evolution at the time-independent frequency uj±. 
Combining expression J3J) with the condition ujq <C 
we find the solution of Eq. (J3J) in the form 



b(t) =l + h(t), 



sin cuot. 



(6) 



The field operator \ can be represented in the usual form 
X = (xo + x') 6xp(— i/ir), xo being the condensate wave 
function and p, being the initial chemical potential. In the 
absence of excitations at T = we can omit the last term 
in Eq. (@J and determine Xo(p) m the universal form. 
Considering excited states, we first linearize Eq. Q in x' ■ 
Under condition luq <C luj_ only long-wave longitudinal 
phonons prove to be involved into the evolution of the 
system. For these phonons, the transverse distribution 
is close to xo (see, e.g., 0). This allows us to simplify 
the equation obtained. Using the second relation in © 
and going over to variable t, we arrive at the following 
equation 



. h d X ' G d* X ' 
dt ~ b 2 (t) [X + X ' 2m Sz 2 



G = U xl (7) 



Involving the symmetry of the problem, we have for the 
operator x' m second quantization 



X 1 = z2xkak; Xk = -j= ■ (t>o(~p), 



(8) 



where &k being the annihilation operator for atoms. At 
luq <C lo±_ the function 4>o is actully close to xo f° r both 
the quasi-lD and Thomas-Fermi cases. In the expression 
(JHJ) we keep the ground state alone for the transverse 
motion assuming that higher states are insignificant in 
the course of evolution. Substitution of the expression 
into Eq.0 gives after the standard averaging over 
the radial variables 



\ + ~2~mJ k + ^"- fe - vGsinLU t(ak + at k ) 



dt \ 

(9) 

G = J d 2 p4>^( p )G( p ). Let us rewrite this equation in 
terms of phonon operators using the usual transforma- 
tion dk = Ukbk + Vkb_ k . When the last term in Eq. © 

is omitted, a set of equations for b k ,b_ k determines the 
well-known Bogoliubov spectrum lo^ (with the replace- 
ment nllo — > G). Considering the general nonstationary 

case, we introduce the substitution bj~ — b k exp(—iu>ot/2) 
and take into account that long times oj$t 3> 1 are most 
interesting for the analysis. Then, within the resonance 
approximation the equations for bk,b_ k are reduced to 
the form 



dbk t A±;i7 i+ 
i—jT = ikb + iE ok bZ k ] 



db + h 

-i-^ = £,kbt k ~ iE kb, (10) 



Eok = (rjG/2h)l(u k + v k )/u k ], (,k = uj k -uj^/2. In fact, we 
consider the evolution of excitations at the background of 
the coherently oscillating condensate. These excitations 
can be found as oscillations of classical Bose-field of the 
condensate (see, e.g., 6j). With this advantage the oper- 
ators in Eqs. (|10fl can be replaced with the classical func- 
tions. If the solution of Eqs. (|10|) is represented in the 
form& fc ,6+ fc ~ exp(aokt), one finds a ok = [\E k\ 2 -t,l] 1/2 ■ 
This result demonstrates the appearance of PR with an 
exponential growth of the phonon occupation numbers 
at l-Eol > |£/c|j which is induced by the modulation of 
the transverse trap frequency. Thus the parametric res- 
onance occurs within the narrow range near lvq with the 
width 2Eq. For sound phonons \(uk +Vk)/uk\ ~ hw k /p 
and for the parameter Eq, we find 
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(11) 



So far we have disregarded the phonon-phonon inter- 
action. With an exponential growth of the phonon num- 
ber the interaction begins to play an essential role. The 
weakness of the phonon-phonon interaction as itself (see 
below) allows us to take only three- and four-phonon 
processes into account. Assuming generation of sound 
phonons alone in the parametric interval 2E , we can 
use the expressions for and obtained within the 
hydrodynamic approximation, see 0. The direct analy- 
sis shows that the contribution of is small in com- 
parison with the term calculated in second order of the 
perturbation theory with respect to owing to small- 
ness of the gas parameter. The dominant term in the 
Hamiltonian of thrcc-phonon interactions has the form 



# (3) = — 



d xvdnv. 



[12) 



Here v, Sn are the operators of velocity and alternating 
part of the density, respectively. In fact, we consider 
the interaction of sound phonons that reduces to 
the one-dimensional problem. Since each three-phonon 
vertex implies the momentum conservation law, at least 
one of three phonons lies beyond the parametric interval 
and, therefore, has zero occupation number at T = 0. 
This makes possible to reduce the expression obtained in 
second order in to the effective Hamiltonian for the 
four-phonon interaction (A(fc) is the Kroneker symbol) 



H, 



HA 



&fc 1 fefe 2 ^3^4A(fc 1 +fc 2 -fc 3 -fc 4 ). (13) 

k\ , . . . ,/C4 



Here all states lie within the interval of 2Eq. Using the 
known expressions for operators v, 6n 7] and condition 
£o <C wo, we obtain A = huj^/pN where N is the total 
number of particles. In addition, second order in H^ 1 
contains the imaginary part related to real decay pro- 
cesses of phonons, which determine the phonon damping. 
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Later, we take into account phenomenologically introduc- 
ing a decrement 7^ . After work one can conclude that 
the parameter j k is small at T = for the geometry un- 
der consideration. Therefore, it is rather easy to satisfy 
the conditions when ctQ k > Ik and the parametric growth 
of phonon number remains. Now let us write the equa- 
tion for b k taking -ff e // into account. Using substitution 

bk = bk exp(—iu>ot/2) and going over to the classical Bosc 
field for phonons, we have 



+ A Kj>k 3 bkA{k + k 2 - h - k A ) 

kik^ki 



(14) 



By means of the equations for bk and b*_ k one can directly 
obtain the equations for the correlators Nk —< b k b k > 
and fk —< bkb-k >• Using conditions 4 « £ « ujq, 
we decouple the arising four-phonon terms within the 
mean- field approximation. As a result, we arrive at a set 
of nonlinear equations that describes the self-consistent 
evolution of interacting phonons in the vicinity of PR 



dt 



=^Nk + E (fk + f*k) + iA(V*fk~rr k ); 



*^=2(-*7+&) fk + 2iE Nk + 2ATN k +4AQfk. 



(15) 



Here Q = Yl ^fc' ! ^ = EA'- Hereafter the evident 

k' v 
relations uj k = ^>-k, Nk = N- k are used. 

In the general case, Eqs. (|15[1 are a set of nonlinear 
integral equations. To simplify the situation, let us con- 
sider the case of finite longitudinal size L, when only a 
single (two-fold degenerate) level lies within the energy 
interval of about 2Eq . The level position £ with regard 
to the resonance center is determined by the relation 
uj k — loq/2 = £. Let us represent the function f k (t) in 
the form f k (t) = \f k (t)\exp[i(f k (t)] = N k (t) exp[iip k (t)} 
(N k > 1). Then, from Eqs. (JT5J we find 



dNk 
dt 

d(fik 
dt 



-2jN k + 2E N k cos ip k ; 



(16) 



-2[£ + AN k + E sin tp k ], A = 3A. 
The stationary solution of this system has the form 



\A\ 



s [l[i -<Pk = T 



E 



(17) 



The signs =p correspond to A ^ 0, respectively. For def- 
initeness, we suppose Eq > 0. The result l|17l) has an 
interesting physical origin. In fact, the interaction -ff e // 
determines the effective renormalization of the phonon 
level, leading to Sujk — AN k for the case concerned. 
Accordingly, £ -»■ £ = £ + AN k , Fro m Eqs. (|T7|) it 
straightforwardly follows that £ s = ± \/ Eq — j 2 . So, at 



Nk = N k the renormalized level reaches the left or right 
edge of the parametric energy interval (within the accu- 
racy of the shift due to 7). As a result, the parametric 
increase of the phonon occupation numbers stops and the 
phase acquires the constant value. The maximal value of 
the phonon number equals 



kmax \A\ 



7,£< E . 



(18) 



The phase is cp k « T^/^ f° r A ^ 0. The phase ip k 
corresponds to the phase correlation of phonon pairs with 
the opposite momenta. The appearance of the phase ipk 
and anomalous correlator fk is connected, first of all, 
with the creation of phonon pairs in the course of the 
evolution as a result of the parametric resonance. The 
phonon-phonon interaction does not disturb the phase 
correlation, and the evolution finishes in the stationary 
condensate of phonon pairs with zero momentum and the 
common phase ip^. 

The character of the temporal evolution, tending 
asymptotically to the values fT7| . l(T%)) . depends essen- 
tially on the relations between the parameters in Eqs. 
(fit))) . First, let 7 = 0. In this case, Eqs. flfi|l have the 
conserved integral of motion 



H = 2E N k smcpk 



AN* 



2£Nk 



(19) 



At the initial time moment when £ = £ + ANk(0) and 
"0 = \JEq— £q, Eas. ((Trj|) imply that cosi^O) = a /E . 
Bearing in mind that sinip k (0) = —^q/Eq, from the sec- 
ond equation in (fit))) we obtain (dtpk/dt)(0) — 0. As 
a consequence, the initial (and conserved) value of Hq 
is equal to Hq w —AN 2 (0) <C Eq- For an arbitrary 
time, Eq. (|19fl can be rewritten as 2Eq sin (p k (t) = 
-[2£ + AN k (t)] - H Q /N k (t). Substituting this relation 
into equations (|l(j[) . one can find 



dN k 
dt 



±2N k a: 



d(fk 
dt 



-ANk 



a = [El - {ANk/2 + £ - H /2Nk) 2 Y/ 2 . The signs ± 
correspond to the regions with |<^fc| < tt/2 and \(fk\ > 
7r/2, respectively. For the region with \(pk\ < 7r/2, the 
solution can be found straightforwardly 



N„(t) 



2t; 



dx 



Af fc (0) 



x[E 2 -(Ax/2 + 2 } 1/2 ' 



(21) 



Here we omitted the small term with |i?o| *C Eq. The 
upper limit of the integral is equal to the value TV™ « 
2(Eq £)/\A\ at which the denominator vanishes. The 
divergence at N k — > is a typical manifestation of the 
parametric resonance which requires the finiteness of the 
initial field amplitude. This can be achieved by taking 
zero-point oscillations into account 0. Owing to these 
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oscillations, we can put N k (0) ~ 1. The time necessary 
for the system to achieve the maximal value N™ is equal 
to 



tr. 



1 

2a~o 



■In 



A 



El 



(22) 



We see that the argument of logarithm is much greater 
than unity at £ < Eq and Eq 3> A. This implies that 
t m l/2ao where l/2ao is the characteristic time of the 
parametric resonance. At t = i m we have |<pfc| = 7r/2. As 
it follows from the second equation in (120ft . at i > t m the 
phase proves to be in the region with \ fk\ > 7r/2. As 
a result, in the first equation the sign becomes negative 
and the phonon number reduces. This is the start of an 
oscillating behaviour. The joint solution of Eqs. (f 2()f) 
has the form of high anharmonic oscillations of both the 
phonon number and the phase around their stationary 
values. At £ = the solution can be written in terms 
of Jacobian elliptic functions. Averaging the dependence 
N k (t) over the large temporal interval t 3> t m , we arrive 
at value AT|. 

The character of evolution changes drastically at 7 7^ 0. 
Namely, a self-averaging arises, leading asymptotically 
to the stationary value iV|. In this case, the integral 
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FIG. 1: ! The phonon numbers Nk(s) versus s — 2~/t. The 
cases (a), (b), and (c) for (Eo/j) = 5, 2, 1.1. The case (d) 
plots Nk(t) versus y(t) — — y>k(t) for (Eo/y) = 5. 

of motion is absent and the solution of Eas. (|16|) should 
be found directly. The damping of the oscillations has 
the decrement close to 7 when ao 7 and t > t m . As 
7 increases, the arrival time to the stationary state de- 
creases as I/7. When 7 is close to ao, the oscillations 
must disappear completely In fact, the time necessary 
for the phonon number to reach the stationary state will 
be determined by Eq. (f 2 21) in which ao is replaced by 
ax = a o — 7- The direct numerical simulation of the sys- 
tem ifTfif demonstrates the picture described. In Figs.l 
(a)-(c), the phonon number Nk as a function of s = 2jt 
is shown at various values of the ratio Eq /j for the fixed 



value of the ratio (Eq/A) = 10 2 and £ = 0. As an illustra- 
tion, Figure 1 (d) displays the phase portrait in the plane 
[Nk, <fk] for (E /j) = 5. Thus the self-consistent tempo- 
ral evolution of interacting phonons near the parametric 
resonance results in the formation of the stationary state 
with the macroscopic phonon number (of the order of 
Eq/A) in a single quantum state. The stationary state 
has one of two energies e s = H(u>q/2± Eq). In the spatial 
sense this is a superposition of two states with the wave 
vectors k s and — k s . 

The involvement of a coherent combination of phonons 
with the momenta fc s and — k s induces the station- 
ary space modulation of the condensate. It is easy 
to reveal determining the particle density as An =< 
X + (p, z)x'{p, z) >, where %' has the form JSJ). Rewriting 
x' in terms of phonon operators and integrating over d 2 p, 
we arrive at the expression for the modulation of the ID 
atomic density 



Sn{z)_2N l 



N 



k ( u l s + v l s ) cos2k s z ~— •^-•cos2fc s z, (23) 



N 



where rS 1 ^ = N/L. We suppose that \5n\/n^ <C 1. 

Let us make here some estimations. As an illustration, 
we consider a quasi-lD gas of sodium atoms. The quasi- 
1D situation occurs under condition p < cjx- Here p ~ 
(an^)Huj± (see 0), a being the 3D scattering length. 



For both = lO^m" 1 

anW « 0.2 and N w 10 4 



and L w 10 cm, we have 
Assuming (luj_/luq) = 5, we 



find 



Nt 



10 4 



\SUJ± 



M 

nW 



3\Sujj_\ 



The above analysis is performed for the case when a sin- 
gle level lies within the energy interval 2Eo. This as- 
sumption imposes certain restrictions on the ratio Eq / Auj 
where Aw = ttc/L is the spacing between the neighbor 
longitudinal levels, c being the sound velocity. Com- 
bining definition of Eq (|llfl and value u>± = 10 1/sec, 
we find (Eq/Aoj) w 10|<5tL>j_|/a;_i_). These estimations 
demonstrate that at |5w^|/cjjJ ~ 10 -2 the conditions 
\5n\/n^ < 1, Eq/Auj < 1 are satisfied and value 
proves to be about 10 2 . In reality, we may increase the 
ratio (\du)±\/uj±). As a consequence, the value iV^ should 
be within interval 10 2 — 10 3 . The parameter j/Eq es- 
sential for the temporal evolution may be changed vary- 
ing the ratio |<5cj^|/ciJj^. Thus, we can conclude that the 
phonon condensation can be realized for a quite realis- 
tic set of parameters. It should be noted that, for the 
finite but sufficiently low temperatures when the initial 
number of phonons with io% ~ ujq is limited, the results 
obtained remain valid. 
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